Chiral spin pairing in helical magnets 
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A concept of chiral spin pairing is introduced to describe a vector-chiral liquid-crystal order in 
frustrated spin systems. It is found that the chiral spin pairing is induced by the coupling to 
phonons through the Dzyaloshinskii-Moriya interaction and the four-spin exchange interaction of 
the Coulomb origin under the edge-sharing network of magnetic and ligand ions. This produces 
two successive second-order phase transitions upon cooling: an 0(2) chiral spin nematic, i.e., spin 
cholesteric, order appears with an either parity, and then the 0(2) symmetry is broken to yield a 
helical magnetic order. Possible candidate materials are also discussed as new multiferroic systems. 

PACS numbers: 75.10.Hk, 05.50.+q, 63.70.+h, 77.80.-c 



The chirality in the electronic spin texture introduced 
by a geometrical frustration and/or the relativistic spin- 
orbit interaction has been one of the key concepts in 
strongly correlated electron systems. The scalar spin chi- 
rality [HQ, i.e., the scalar triple product of three non- 
coplanar spins, is odd under the time-reversal (T) and 
even under the space-inversion (I), and gives rise to a 
large anomalous Hall effect @, S IE 0, Q ■ On the other 
hand, central to this Letter is the vector spin chirality Q , 
i.e., the vector product of two noncollinear spins. It is 
T-even and /-odd, and can produce the ferroelectric po- 
larization in Mott insulators through the spin-orbit in- 
teraction 
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The vector-chiral spin order is realized in conventional 
helical magnets. In principle, it is even possible that the 
chiral or parity symmetry is broken but the time-reversal 
symmetry is not. This state having hnite spin correlation 
lengths is categorized into a liquid or a liquid crystal of 
spins, which is of our main interest, while the magneti- 
cally ordered state into a solid. The issue of the chiral 
spin order in the absence of any magnetic order is not 
only a problem of statistical mechanics, but has been in- 
tensively studied since the discovery of high-T c cuprates, 
and even gives a new "multiferroic" phenomenon which 
has attracted current great interests, as discussed below. 
However, conditions for the order being stabilized have 
not been fully understood. 

It has been ar gued that the vector or pseudoscalar chi- 
ral ordering 21, 22] takes place separately from the 
XY transition in the antiferromagnetic (AF) XY model 
on the triangular lattice [23j |. The intermediate phase is 
characterized by the vector-chiral order. It has also been 
obtained for frustrated quantum spin systems in one di- 
mension [lij], and discussed in two dimensions [25[ and 
for odd-time states [26| . However, it usually appears only 
in a tiny region of the global phase diagram. For classical 
AF Heisenberg and XY models on stacked triangular lat- 
tices, Monte-Carlo simulations and field-theoretical anal- 
yses have suggested a single phase transition from param- 
agnet to helical magnet, which is either weakly first-order 



or a second-order one, which belongs to a different uni- 
versality class from the O(N) model [27|, EE SU • Then, 
the chiral spin liquid-crystal phase does not appear. 

In this Letter, we introduce a picture of chiral spin 
pairing for the vector-chiral spin order which has an anal- 
ogy to a cholesteric liquid crystal defined as a chiral ne- 
matic [30j, and propose a new strategy to realize this 
phase. It is found that when the coupling of spins to opti- 
cal phonons through the Dzyaloshinskii-Moriya (DM) in- 
teraction and/or the Coulombic four-spin ring-exchange 
interaction under the edge-sharing network of magentic 
and ligand ions are robust against the magnetostriction, 
these interactions bind spin pairs with the vector chiral- 
ity without any magnetic long-range order. Then, upon 
cooling, the transition from paramagnet to helical mag- 
net takes place through two steps of lowering the sym- 
metry. Namely, an 0(2) spin cholesteric and a helical 
magnetic phases succesively appear. (See Fig. [I]). 

We begin with the Ginzbrug-Landau Hamiltonian for 
incommensurate helical Heisenberg magnets with the or- 
dering wave vector Q, which can be expressed in terms 
of two independent long wavelength modes of the 0(3) 
field S r for the spin at a site r, namely, S q = S r e~ lq r 
with q ~ ±Q, or equivalently the long-wavelength modes 

Sq = ^{^Q+q + S-Q+q) an d b q = —^(Sq+ q — S-Q +q ) 

with |q| < A and the momentum cutoff A. This gives 



S r w a r cos Q ■ r b r sin Q ■ r, 



(1) 



in terms of the inverse Fourier transforms (a r ,b r ) — 
a o J A (2Ti) d ("<J' bq)e lq r with the lattice constant ao- Here, 
"ps" represents a weak equality except the irrelevant 
rapidly varying fields with |q| > A. Substituting Eq. (TT]) 
into the Heisenberg Hamiltonian, one obtains the contin- 
uum model [27| given by 
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We have neglected the spatial anisotropy in the gradient 
terms for simplicity, and rescaled the variables so that 



2 



coefficients to those terms become unity. Quartic inter- 
action terms with the coupling constants u — uq > 
and v — Vq — —Uq/3 can be obtained by softening 
the constraint that amplitudes of spins must be fixed, 
S(r) 2 w S 2 , and substituting Eq. (JTJ) into the quar- 
tic term ^ Jdr(S 2 ) 2 . Then, the Hamiltonian given by 
Eq. |(3J) is isomorphic to a classical 0(3) x 0(2) nonlincar- 
a model [27[. Now helical and collinear magnetic orders 
are expressed by (a) x (b) ^ and (a) x (b) = 0, re- 
spectively. The vector-chiral order is characterized by 
(a x b) ^ 0. u and v are modified from the Heisenberg 
values by additional interactions as shown later. 
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FIG. 1: (Color online) Phase diagram of the Ginzburg- 
Landau Hamiltonian given by Eq. ^ and the ordering pat- 
terns. Note that ft, = — 8u. For details, see the text. 

Figure [1] presents the main result of this paper, i.e., the 
phase diagram for d = 3 obtained in the space of (ft, v/u) 
for the model Eq. ©. Here we introduce dimensionless 
coupling constants u = uciQ/(2ir 2 A) and v = va§/(27r 2 A), 
as well as ft = p/A 2 which corresponds to the tempera- 
ture. The spin cholesteric, i.e., chiral spin pairing state, 
appears above the helical magnetic ordering transition 
temperature when v/u is less than -1/3. When ft is 
negatively large, a negative/positive v stabilizes a heli- 
cal/collinear magnetic order. These two phases are sep- 
arated by the first-order phase transition at v = 0. 

To derive this result, we have employed the mode cou- 
pling approximation, where the average value and the op- 
timal Gaussian fluctuation form of the order parameters 



are determined by the variational principle [34| . First, 
we decompose the fields into the condensed and normal 
components as (a q , b q ) = S q (A, B) + (5a q , 6b q ). (a r ) = A 
and (b r ) = B give the ordered mangetic moment through 
Eq. (fT]). (A,B) is subject to the stationary condition 
of the free energy with respect to (Set, Sb) . To treat 
the Gaussian fluctuations around the saddle point, we 
consider the variational quadratic Hamiltonian TL var = 



(ao/2n)- 3 JdqJ2i 1=1 (Sa q ,Sb q ) 



q 



'(SaLvSbLa) 



with the 6 x 6-matrix Green's function G q defined by 



(S hJ q 2 +^)p°+S^p z 



•A^ + AlV (3) 



with the 2x2 identity matrix p° and the Pauli matrices 
(p x ,p y , p z ) operating on the (a, b) space. Here, the varia- 
tional parameters , Sp^ , A l J , and A l J are determined 
from the self-consistent equations, 
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(u + 3v) (A l B j - A>B l + irf) , 



(4a) 
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where 7 ^ = (Sa\Sa 3 r + Sb\5b> r ) , 7 * : 
v« = (6a*.8a£ - Sb\.5V r ), T « 

v % 3 = 
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i"f l J ee (5a\.5b> r — Sb l r Sal) are calculated from 



Sbi + SbiSai), 
- SbtSai), and 



in 
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(2^y 



(5) 



where the trace Tr is taken only in the (a, b) space. 

We start from the disordered phase where all the 
variational parameters including (A, B) vanish except 
the mass p r = p m . Then, the self-consistent equa- 
tions (U)) can be linearized with respect to the dimension- 
less variational order parameters /2 y = /U y /A 2 (i ^ j), 
(Sfti lj , A\ 3 ) ee (5^,Af)/A 2 , and Ajf = A^/A 2 ; 



ft'' 
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p r / A 2 = ft + 4(2u + v)g(ft r ), (6a) 

-2(u-v)ft^f(ft r ) for i^j, (6b) 

-2 ( %j ) [-3*% + (u + v)} f(ft r ), (6c) 

-2A^(u + 3v)f(p r ), (6d) 



with 



g(ft r ) = 1 — y/i r arctan 



1 1 1 

arctan ■ 



ft r 1 + ft r 



(7a) 
(7b) 



If we ignore the spin pairing, i.e., <5/i y = A Z J = A l J = 0, 
Eq. (|6aj) always gives a direct phase transition from the 
paramagnet to a helical (collinear) magnet for v < 
(v > 0) at a critical point ft r = 0. However, a serious 
consideration on the spin pairs Sp^ , A% , and A l J reveals 
that a spin liquid-crystal phase emerges before the mag- 
netic order sets in. Note that g(ft r ) and f(ft r ) are positive 
and that f(p r ) diverges as 7r/(4v / /S F ) for fti r — > 0. Then, 
as far as the Hamiltonian is stable, i.e., u + v > 0, there 
appear two possibilities for spin paired states, (i) When 
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the interaction for the antisymmetric spin pairing chan- 
nel is attractive, namely, v/u = v/u < —1/3 in Eq. (|6d[) . 
then it forms the chiral spin pairs A* J ' ^ through a 
second-order phase transition. The critical value of u r is 
determined by /(/*„) = — l/2(u+35) from Eq. (|6d|) with 
/z£ r > 0: the magnetic long-range order is absent. This 
represents a 0(2)-symmetric chiral spin nematic (26j or 
spin cholesteric. (ii) When v/u = v/u > 1/2 in Eq. ([6c]) . 
the symmetric spin pairs are formed, (Su ll ,Al l ) 7^ 0. 
However, this condition is usually difficult to be realized 
starting from the Heisenberg point v/u = — 1/3. In the 
following, we concentrate on the case of (i) v/u < —1/3, 
which is more likely to occur. 

In the spin cholesteric phase characterized by 
^■a Hy e yz with the fully antisymmetric tensor e+ji, 
the Green's function can be described by 
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-q 2 -a{l-S lz )A a 



(8) 



Here, = ^ xa: = = ^ zz , and A a are determined 

by the numerical solution to the self-consistent equations, 

/x_L = fl + 2uh+(fl ± ,0) + 2{u + i!)h + (fl h A a ), (9a) 
A|| = A+(« + «)^ + (A±,0) + (3u + 5)/i + (/i||,A a ),(9b) 
A a - -(u + 3«)/i-(/i||,A ), (9c) 



which are derived from Eqs. (|4a|) and (|4dj) . where 
^(/i, A a ) = A a )±p(/I + A ). Since || < fix holds 
generally, the spins that are originally 0(3)-symmetric 
now obtain the easy-plane XY anisotropy. Accordingly, 
the spin correlation has three nonvanishing branches; 



(StSZ) = Xr,r'(/U) cos Q-(r-r'), 

(S r S r i) = — / A.r,r' VS|| ff 
cr=± 



(10) 



^ Xrv'(^ 2 ) cos Q-(r-r'), (11) 



with i = x,y and the chiral spin correlation 

; ll 



(S?S£-£*££) = J] ^,r'(C) sinQ-(r-r'), (12) 



CT=zt 



where £ CT = 1 / — oA a and 



a 3 A 



"^"Tl 77 ex P ( _ £l r - r 'l) 

47r 2 A|r-r'| v 1 



(13) 



for a long-distance decay. These spin correlations are all 
short-ranged with two correlation lengths ~ £± . The 
sign of the chirality A a can be determined from that of 
(S^S^, —S%S*/)- Further decreasing u, fin becomes equal 
to A a , leading to a divergence of either £ + or and thus 
the second-order phase transition to the helical magnet. 
Note that all the above structures should appear in the 
static spin correlations, which can be directly studied 
by the polarization dependence of the incident and/or 
scattered neutrons [35j |. 



In fact, the coupling constants u and v are modified by 
(i) the coupling of spins to longitudinal phonons with the 
wavevector ±2Q through the magnetostriction, (ii) that 
to transverse phonons with the wavevector through 
the DM interaction, and (iii) the four-spin ring-exchange 
interactions, i.e., u = uq + 5u ms + 5u + <5u rms and 
v = -uo/3 + Sv ms + Sv DM + <5v rillg . 

(i) The magnetostriction arises from 



n . n . 



p 2 



1 



2M ' 2 KX * 



(14) 

with the spatial derivative 5 J of the exchange interaction, 
the unit vector e r t- r directing from r to r', the mass 
M of the magnetic ion, the spring constant K associated 
with its shift X ri and the conjugate momentum P r . Sub- 
stituting Eq. |T]) into Eq. (fT4|) and integrating over X r , 
we obtain 5u MS = 5v MS = (cr 2 SJsinQ ■ S) 2 /AK where 
c is the average spin-wave velocity. 

(ii) Additional terms to the Hamiltonian reads 
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-D r 
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Vt.t' 

2m l 



1 ^ 2 
~h K,X T . 



(15) 



Here, D r . 



\x T 



X 6f 



represents the anti- 



symmetric DM coupling vector 31| between the nearest- 
neighbor spin pair S r and S r < located at the discrete 
lattice sites r and r', respectively, with the coupling con- 
stant A and a transverse shift x r . r ' of the ligand ion [l7| . 
For simplicity, we have assumed that the ligand ion is 
located at the center of nearest-neighbor bonds of spins. 
Because of this coupling, the chiral order (a X b) 7^ 
simultaneously produces both macroscopic averages of x 
and the electric polarization in proportion to (a x b) and 
Q x ((ox b)), respectively. Note that applied electric field 
plays the same role as the average value of x. Substitut- 
ing Eq. |T]) into Eq. (fl"5|) and integrating over x r y, we 
obtain 5u DM = and Sv DM = -(c~ 2 A sin Q ■ S) 2 /2k with 
the nearest-neighbor coordination vector S. 

(iii) In the case of corner-sharing network of mag- 
netic and ligand ions, four-spin ring-exchange interaction 
arises from the strong-coupling expansion. Substituting 
Eq. |l| into the ring-exchange Hamiltonian aorund four 
sites ri (i = 1, • • • , 4) in a xy-plane plaquette, 



J ring (P ru 



P 



V2 ,t*3 ;r*4,ri 



p 



ri ,r 3 ;r 2 ,7" 4 ; 



(16) 



with P r 



{S r ■ S r >)(Sr" • S r ">) for spin-1/2 [32, 



M, we obtain <5u rmg 



1 J rmg (cos 2 Q x + cos 2 Q y - 1/2) 



and 6v Tlne = c" 4 J nng (sin 2 Q x + sin 2 Q y - 1/2). 

On the other hand, in the case of edge-sharing net- 
work, the above ring-exchange is suppressed, while the 
second-order strong-coupling perturbation of the four- 
site Coulomb interaction in the unit square gives a still 
fnitc contribution of the form 



-^1 [-^ri ,V2 ;t 



Pt 



K,P, 



ri ,i" 3 ;r 2 ,r 4 



(17) 



4 



4(W - /i 243 ) 2 /^ < Ki, 
>(P - n)(j){p' - rj)(j){p' - 



where K x ~ 4I 2 234 /t/, K 2 

and = J dpdp' \^L p i | v / vp — 'iivkh ~ 'j; 

rk)4>{p — re) with the local Coulomb repulsion t/ and 
the Wannier function <f>{p) centered at the magnetic ion 
site. Again substituting Eq. fl} into Eq. (fT"7|) . we obtain 
^■ing = -0-4(2^ + X 2 )(cos 2 Q x + cos 2 Q y - 1/2) and 

^ring = c -4 Xl ( cos 2 + CQS 2 Q y -\) + K 2 /2. 

Putting the above considerations together, the devi- 
ation of the parameter v/u from the Heisenberg value 
— 1/3, which plays a crucial role in determining the phase 
diagram, is eventually given by 
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y-v sin Q ■ 8 
s 



3c 4 uq 



2K 



3A 2 . 
2k 



Namely, the DM interaction and the Coulombic four-spin 
ring-exchange yield a negative shift of u/v and favor the 
chiral spin pairing, whereas the magnetostriction and the 
kinetic four-spin ring-exchange interaction do not. 

Let us discuss a possibility of observing the spin 
cholesteric phase in real systems. First, the easy-axis 
single spin anisotropy is unfavorable, since it produces a 
mass between the two fields to be paired and suppresses 
the pairing. On the other hand, the chiral order remains 
in the easy-plane anisotropy since the chiral order pa- 
rameter is now a scalar which forms a long-range order 
more easily than the 0(3) vector. In the case of i?Mn03 
(R=Tb, Dy)[l3], the spin of the Mn ions is S = 2 and 
the single spin anisotropy plays a dominant role for the 
successive phase transitions from paramagnet to incom- 
mensurate collinear states and to incommensurate helical 
state as the temperature is lowered. Especially the tran- 
sition to the helical state is of first-order, and hence the 
spin-pairing above the transition is unlikely. Therefore 
S = 1/2 is preferable. Then, the single spin anisotropy 
is forbidden, while the anisotropy of the exchange in- 
teraction is still possible. The edge-sharing network of 
magnetic and ligand ions is also advantageous over the 
corner-sharing, since it suppresses the magnetostriction 
and the kinetic four-spin ring-exchange interaction. Cu 
ion with spin S = 1/2 is an interesting candidate from 
this viewpoint. Actually multiferroic behaviors have been 
recently discovered in edge-sharing LiCuV02 [H| and 
LiCu202 [HI]. Interestingly, the phase transitions to the 
helical state in both materials seem to be second order. 
A recent ESR experiment shows that the DM interaction 
of the order of 1 K survives [36j] , It is highly desirable 
to look for this intriguing possibility of chiral spin pair- 
ing by careful polarized neutron-scattering and ferroelec- 
tric polarization measurements slightly above the helical 
magnetic ordering temperature. 
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